Vitamin D regulates microbiome-dependent cancer immunity.

A role for vitamin D in immune modulation and in cancer has been suggested. In this work, we report that mice with increased availability of vitamin D display greater immune-dependent resistance to transplantable cancers and augmented responses to checkpoint blockade immunotherapies. Similarly, in humans, vitamin D-induced genes correlate with improved responses to immune checkpoint inhibitor treatment as well as with immunity to cancer and increased overall survival. In mice, resistance is attributable to the activity of vitamin D on intestinal epithelial cells, which alters microbiome composition in favor of Bacteroides fragilis, which positively regulates cancer immunity. Our findings indicate a previously unappreciated connection between vitamin D, microbial commensal communities, and immune responses to cancer. Collectively, they highlight vitamin D levels as a potential determinant of cancer immunity and immunotherapy success.

High-throughput anaerobic screening for identifying compounds acting against gut bacteria in monocultures or communities.

The human gut microbiome is a key contributor to health, and its perturbations are linked to many diseases. Small-molecule xenobiotics such as drugs, chemical pollutants and food additives can alter the microbiota composition and are now recognized as one of the main factors underlying microbiome diversity. Mapping the effects of such compounds on the gut microbiome is challenging because of the complexity of the community, anaerobic growth requirements of individual species and the large number of interactions that need to be quantitatively assessed. High-throughput screening setups offer a promising solution for probing the direct inhibitory effects of hundreds of xenobiotics on tens of anaerobic gut bacteria. When automated, such assays enable the cost-effective investigation of a wide range of compound-microbe combinations. We have developed an experimental setup and protocol that enables testing of up to 5,000 compounds on a target gut species under strict anaerobic conditions within 5 d. In addition, with minor modifications to the protocol, drug effects can be tested on microbial communities either assembled from isolates or obtained from stool samples. Experience in working in an anaerobic chamber, especially in performing delicate work with thick chamber gloves, is required for implementing this protocol. We anticipate that this protocol will accelerate the study of interactions between small molecules and the gut microbiome and provide a deeper understanding of this microbial ecosystem, which is intimately intertwined with human health.

Metabolic exchanges are ubiquitous in natural microbial communities.

Microbial communities drive global biogeochemical cycles and shape the health of plants and animals-including humans. Their structure and function are determined by ecological and environmental interactions that govern the assembly, stability and evolution of microbial communities. A widely held view is that antagonistic interactions such as competition predominate in microbial communities and are ecologically more important than synergistic interactions-for example, mutualism or commensalism. Over the past decade, however, a more nuanced picture has emerged, wherein bacteria, archaea and fungi exist within interactive networks in which they exchange essential and non-essential metabolites. These metabolic interactions profoundly impact not only the physiology, ecology and evolution of the strains involved, but are also central to the functioning of many, if not all, microbiomes. Therefore, we advocate for a balanced view of microbiome ecology that encompasses both synergistic and antagonistic interactions as key forces driving the structure and dynamics within microbial communities.

Modelling of microbial interactions in anaerobic digestion: from black to glass box.

Anaerobic and microaerophilic environments are pervasive in nature, providing essential contributions to the maintenance of human health, biogeochemical cycles and the Earth's climate. These ecological niches are characterised by low free oxygen and oxidants, or lack thereof. Under these conditions, interactions between species are essential for supporting the growth of syntrophic species and maintaining thermodynamic feasibility of anaerobic fermentation. Kinetic models provide a simplified view of complex metabolic networks, while genome-scale metabolic models and flux-balance analysis (FBA) aim to unravel these systems as a whole. The target of this review is to outline the main similarities, differences and challenges associated with kinetic and metabolic modelling, and describe state-of-the-art modelling practices for studying syntrophies in the anaerobic digestion (AD) case study.

A new metabolic model of Drosophila melanogaster and the integrative analysis of Parkinson's disease.

High conservation of the disease-associated genes between flies and humans facilitates the common use of Drosophila melanogaster to study metabolic disorders under controlled laboratory conditions. However, metabolic modeling studies are highly limited for this organism. We here report a comprehensively curated genome-scale metabolic network model of Drosophila using an orthology-based approach. The gene coverage and metabolic information of the draft model derived from a reference human model were expanded via Drosophila-specific KEGG and MetaCyc databases, with several curation steps to avoid metabolic redundancy and stoichiometric inconsistency. Furthermore, we performed literature-based curations to improve gene-reaction associations, subcellular metabolite locations, and various metabolic pathways. The performance of the resulting Drosophila model (8,230 reactions, 6,990 metabolites, and 2,388 genes), iDrosophila1 (https://github.com/SysBioGTU/iDrosophila), was assessed using flux balance analysis in comparison with the other currently available fly models leading to superior or comparable results. We also evaluated the transcriptome-based prediction capacity of iDrosophila1, where differential metabolic pathways during Parkinson's disease could be successfully elucidated. Overall, iDrosophila1 is promising to investigate system-level metabolic alterations in response to genetic and environmental perturbations.

Unravelling metabolic cross-feeding in a yeast-bacteria community using 13 C-based proteomics.

Cross-feeding is fundamental to the diversity and function of microbial communities. However, identification of cross-fed metabolites is often challenging due to the universality of metabolic and biosynthetic intermediates. Here, we use 13 C isotope tracing in peptides to elucidate cross-fed metabolites in co-cultures of Saccharomyces cerevisiae and Lactococcus lactis. The community was grown on lactose as the main carbon source with either glucose or galactose fraction of the molecule labelled with 13 C. Data analysis allowing for the possible mass-shifts yielded hundreds of peptides for which we could assign both species identity and labelling degree. The labelling pattern showed that the yeast utilized galactose and, to a lesser extent, lactic acid shared by L. lactis as carbon sources. While the yeast provided essential amino acids to the bacterium as expected, the data also uncovered a complex pattern of amino acid exchange. The identity of the cross-fed metabolites was further supported by metabolite labelling in the co-culture supernatant, and by diminished fitness of a galactose-negative yeast mutant in the community. Together, our results demonstrate the utility of 13 C-based proteomics for uncovering microbial interactions.

Operation of a TCA cycle subnetwork in the mammalian nucleus.

Nucleic acid and histone modifications critically depend on the tricarboxylic acid (TCA) cycle for substrates and cofactors. Although a few TCA cycle enzymes have been reported in the nucleus, the corresponding pathways are considered to operate in mitochondria. Here, we show that a part of the TCA cycle is operational also in the nucleus. Using 13C-tracer analysis, we identified activity of glutamine-to-fumarate, citrate-to-succinate, and glutamine-to-aspartate routes in the nuclei of HeLa cells. Proximity labeling mass spectrometry revealed a spatial vicinity of the involved enzymes with core nuclear proteins. We further show nuclear localization of aconitase 2 and 2-oxoglutarate dehydrogenase in mouse embryonic stem cells. Nuclear localization of the latter enzyme, which produces succinyl-CoA, changed from pluripotency to a differentiated state with accompanying changes in the nuclear protein succinylation. Together, our results demonstrate operation of an extended metabolic pathway in the nucleus, warranting a revision of the canonical view on metabolic compartmentalization.

Bioactivity assessment of natural compounds using machine learning models trained on target similarity between drugs.

Natural compounds constitute a rich resource of potential small molecule therapeutics. While experimental access to this resource is limited due to its vast diversity and difficulties in systematic purification, computational assessment of structural similarity with known therapeutic molecules offers a scalable approach. Here, we assessed functional similarity between natural compounds and approved drugs by combining multiple chemical similarity metrics and physicochemical properties using a machine-learning approach. We computed pairwise similarities between 1410 drugs for training classification models and used the drugs shared protein targets as class labels. The best performing models were random forest which gave an average area under the ROC of 0.9, Matthews correlation coefficient of 0.35, and F1 score of 0.33, suggesting that it captured the structure-activity relation well. The models were then used to predict protein targets of circa 11k natural compounds by comparing them with the drugs. This revealed therapeutic potential of several natural compounds, including those with support from previously published sources as well as those hitherto unexplored. We experimentally validated one of the predicted pair's activities, viz., Cox-1 inhibition by 5-methoxysalicylic acid, a molecule commonly found in tea, herbs and spices. In contrast, another natural compound, 4-isopropylbenzoic acid, with the highest similarity score when considering most weighted similarity metric but not picked by our models, did not inhibit Cox-1. Our results demonstrate the utility of a machine-learning approach combining multiple chemical features for uncovering protein binding potential of natural compounds.

Molecular biology for green recovery-A call for action.

Molecular biology holds a vast potential for tackling climate change and biodiversity loss. Yet, it is largely absent from the current strategies. We call for a community-wide action to bring molecular biology to the forefront of climate change solutions.

Co-translational assembly orchestrates competing biogenesis pathways.

During the co-translational assembly of protein complexes, a fully synthesized subunit engages with the nascent chain of a newly synthesized interaction partner. Such events are thought to contribute to productive assembly, but their exact physiological relevance remains underexplored. Here, we examine structural motifs contained in nucleoporins for their potential to facilitate co-translational assembly. We experimentally test candidate structural motifs and identify several previously unknown co-translational interactions. We demonstrate by selective ribosome profiling that domain invasion motifs of beta-propellers, coiled-coils, and short linear motifs may act as co-translational assembly domains. Such motifs are often contained in proteins that are members of multiple complexes (moonlighters) and engage with closely related paralogs. Surprisingly, moonlighters and paralogs assemble co-translationally in only some but not all of the relevant biogenesis pathways. Our results highlight the regulatory complexity of assembly pathways.

Unravelling the collateral damage of antibiotics on gut bacteria.

Antibiotics are used to fight pathogens but also target commensal bacteria, disturbing the composition of gut microbiota and causing dysbiosis and disease1. Despite this well-known collateral damage, the activity spectrum of different antibiotic classes on gut bacteria remains poorly characterized. Here we characterize further 144 antibiotics from a previous screen of more than 1,000 drugs on 38 representative human gut microbiome species2. Antibiotic classes exhibited distinct inhibition spectra, including generation dependence for quinolones and phylogeny independence for ß-lactams. Macrolides and tetracyclines, both prototypic bacteriostatic protein synthesis inhibitors, inhibited nearly all commensals tested but also killed several species. Killed bacteria were more readily eliminated from in vitro communities than those inhibited. This species-specific killing activity challenges the long-standing distinction between bactericidal and bacteriostatic antibiotic classes and provides a possible explanation for the strong effect of macrolides on animal3-5 and human6,7 gut microbiomes. To mitigate this collateral damage of macrolides and tetracyclines, we screened for drugs that specifically antagonized the antibiotic activity against abundant Bacteroides species but not against relevant pathogens. Such antidotes selectively protected Bacteroides species from erythromycin treatment in human-stool-derived communities and gnotobiotic mice. These findings illluminate the activity spectra of antibiotics in commensal bacteria and suggest strategies to circumvent their adverse effects on the gut microbiota.

Adaptive laboratory evolution of microbial co-cultures for improved metabolite secretion.

Adaptive laboratory evolution has proven highly effective for obtaining microorganisms with enhanced capabilities. Yet, this method is inherently restricted to the traits that are positively linked to cell fitness, such as nutrient utilization. Here, we introduce coevolution of obligatory mutualistic communities for improving secretion of fitness-costly metabolites through natural selection. In this strategy, metabolic cross-feeding connects secretion of the target metabolite, despite its cost to the secretor, to the survival and proliferation of the entire community. We thus co-evolved wild-type lactic acid bacteria and engineered auxotrophic Saccharomyces cerevisiae in a synthetic growth medium leading to bacterial isolates with enhanced secretion of two B-group vitamins, viz., riboflavin and folate. The increased production was specific to the targeted vitamin, and evident also in milk, a more complex nutrient environment that naturally contains vitamins. Genomic, proteomic and metabolomic analyses of the evolved lactic acid bacteria, in combination with flux balance analysis, showed altered metabolic regulation towards increased supply of the vitamin precursors. Together, our findings demonstrate how microbial metabolism adapts to mutualistic lifestyle through enhanced metabolite exchange.

Towards a mechanistic understanding of reciprocal drug-microbiome interactions.

Broad-spectrum antibiotics target multiple gram-positive and gram-negative bacteria, and can collaterally damage the gut microbiota. Yet, our knowledge of the extent of damage, the antibiotic activity spectra, and the resistance mechanisms of gut microbes is sparse. This limits our ability to mitigate microbiome-facilitated spread of antibiotic resistance. In addition to antibiotics, non-antibiotic drugs affect the human microbiome, as shown by metagenomics as well as in vitro studies. Microbiome-drug interactions are bidirectional, as microbes can also modulate drugs. Chemical modifications of antibiotics mostly function as antimicrobial resistance mechanisms, while metabolism of non-antibiotics can also change the drugs' pharmacodynamic, pharmacokinetic, and toxic properties. Recent studies have started to unravel the extensive capacity of gut microbes to metabolize drugs, the mechanisms, and the relevance of such events for drug treatment. These findings raise the question whether and to which degree these reciprocal drug-microbiome interactions will differ across individuals, and how to take them into account in drug discovery and precision medicine. This review describes recent developments in the field and discusses future study areas that will benefit from systems biology approaches to better understand the mechanistic role of the human gut microbiota in drug actions.

Metabolic cooperation and spatiotemporal niche partitioning in a kefir microbial community.

Microbial communities often undergo intricate compositional changes yet also maintain stable coexistence of diverse species. The mechanisms underlying long-term coexistence remain unclear as system-wide studies have been largely limited to engineered communities, ex situ adapted cultures or synthetic assemblies. Here, we show how kefir, a natural milk-fermenting community of prokaryotes (predominantly lactic and acetic acid bacteria) and yeasts (family Saccharomycetaceae), realizes stable coexistence through spatiotemporal orchestration of species and metabolite dynamics. During milk fermentation, kefir grains (a polysaccharide matrix synthesized by kefir microorganisms) grow in mass but remain unchanged in composition. In contrast, the milk is colonized in a sequential manner in which early members open the niche for the followers by making available metabolites such as amino acids and lactate. Through metabolomics, transcriptomics and large-scale mapping of inter-species interactions, we show how microorganisms poorly suited for milk survive in-and even dominate-the community, through metabolic cooperation and uneven partitioning between grain and milk. Overall, our findings reveal how inter-species interactions partitioned in space and time lead to stable coexistence.

Proteomic characterization of extracellular vesicles produced by several wine yeast species.

In winemaking, the use of alternative yeast starters is becoming increasingly popular. They contribute to the diversity and complexity of wine sensory features and are typically used in combination with Saccharomyces cerevisiae, to ensure complete fermentation. This practice has drawn the interest on interactions between different oenological yeasts, which are also relevant in spontaneous and conventional fermentations, or in the vineyard. Although several interactions have been described and some mechanisms have been suggested, the possible involvement of extracellular vesicles (EVs) has not yet been considered. This work describes the production of EVs by six wine yeast species (S. cerevisiae, Torulaspora delbrueckii, Lachancea thermotolerans, Hanseniaspora uvarum, Candida sake and Metschnikowia pulcherrima) in synthetic grape must. Proteomic analysis of EV-enriched fractions from S. cerevisiae and T. delbrueckii showed enrichment in glycolytic enzymes and cell-wall-related proteins. The most abundant protein found in S. cerevisiae, T. delbrueckii and L. thermotolerans EV-enriched fractions was the enzyme exo-1,3-ß-glucanase. However, this protein was not involved in the here-observed negative impact of T. delbrueckii extracellular fractions on the growth of other yeast species. These findings suggest that EVs may play a role in fungal interactions during wine fermentation and other aspects of wine yeast biology.

Yeast Genome-Scale Metabolic Models for Simulating Genotype-Phenotype Relations.

Understanding genotype-phenotype dependency is a universal aim for all life sciences. While the complete genotype-phenotype relations remain challenging to resolve, metabolic phenotypes are moving within the reach through genome-scale metabolic model simulations. Genome-scale metabolic models are available for commonly investigated yeasts, such as model eukaryote and domesticated fermentation species Saccharomyces cerevisiae, and automatic reconstruction methods facilitate obtaining models for any sequenced species. The models allow for investigating genotype-phenotype relations through simulations simultaneously considering the effects of nutrient availability, and redox and energy homeostasis in cells. Genome-scale models also offer frameworks for omics data integration to help to uncover how the translation of genotypes to the apparent phenotypes is regulated at different levels. In this chapter, we provide an overview of the yeast genome-scale metabolic models and the simulation approaches for using these models to interrogate genotype-phenotype relations. We review the methodological approaches according to the underlying biological reasoning in order to inspire formulating novel questions and applications that the genome-scale metabolic models could contribute to. Finally, we discuss current challenges and opportunities in the genome-scale metabolic model simulations.

Antibiotics-induced monodominance of a novel gut bacterial order.

OBJECTIVE: The composition of the healthy human adult gut microbiome is relatively stable over prolonged periods, and representatives of the most highly abundant and prevalent species have been cultured and described. However, microbial abundances can change on perturbations, such as antibiotics intake, enabling the identification and characterisation of otherwise low abundant species. DESIGN: Analysing gut microbial time-series data, we used shotgun metagenomics to create strain level taxonomic and functional profiles. Community dynamics were modelled postintervention with a focus on conditionally rare taxa and previously unknown bacteria. RESULTS: In response to a commonly prescribed cephalosporin (ceftriaxone), we observe a strong compositional shift in one subject, in which a previously unknown species, UBorkfalki ceftriaxensis, was identified, blooming to 92% relative abundance. The genome assembly reveals that this species (1) belongs to a so far undescribed order of Firmicutes, (2) is ubiquitously present at low abundances in at least one third of adults, (3) is opportunistically growing, being ecologically similar to typical probiotic species and (4) is stably associated to healthy hosts as determined by single nucleotide variation analysis. It was the first coloniser after the antibiotic intervention that led to a long-lasting microbial community shift and likely permanent loss of nine commensals. CONCLUSION: The bloom of UB. ceftriaxensis and a subsequent one of Parabacteroides distasonis demonstrate the existence of monodominance community states in the gut. Our study points to an undiscovered wealth of low abundant but common taxa in the human gut and calls for more highly resolved longitudinal studies, in particular on ecosystem perturbations.

Metabolic models and gene essentiality data reveal essential and conserved metabolism in prokaryotes.

Essential metabolic reactions are shaping constituents of metabolic networks, enabling viable and distinct phenotypes across diverse life forms. Here we analyse and compare modelling predictions of essential metabolic functions with experimental data and thereby identify core metabolic pathways in prokaryotes. Simulations of 15 manually curated genome-scale metabolic models were integrated with 36 large-scale gene essentiality datasets encompassing a wide variety of species of bacteria and archaea. Conservation of metabolic genes was estimated by analysing 79 representative genomes from all the branches of the prokaryotic tree of life. We find that essentiality patterns reflect phylogenetic relations both for modelling and experimental data, which correlate highly at the pathway level. Genes that are essential for several species tend to be highly conserved as opposed to non-essential genes which may be conserved or not. The tRNA-charging module is highlighted as ancestral and with high centrality in the networks, followed closely by cofactor metabolism, pointing to an early information processing system supplied by organic cofactors. The results, which point to model improvements and also indicate faults in the experimental data, should be relevant to the study of centrality in metabolic networks and ancient metabolism but also to metabolic engineering with prokaryotes.

Fast automated reconstruction of genome-scale metabolic models for microbial species and communities.

Genome-scale metabolic models are instrumental in uncovering operating principles of cellular metabolism, for model-guided re-engineering, and unraveling cross-feeding in microbial communities. Yet, the application of genome-scale models, especially to microbial communities, is lagging behind the availability of sequenced genomes. This is largely due to the time-consuming steps of manual curation required to obtain good quality models. Here, we present an automated tool, CarveMe, for reconstruction of species and community level metabolic models. We introduce the concept of a universal model, which is manually curated and simulation ready. Starting with this universal model and annotated genome sequences, CarveMe uses a top-down approach to build single-species and community models in a fast and scalable manner. We show that CarveMe models perform closely to manually curated models in reproducing experimental phenotypes (substrate utilization and gene essentiality). Additionally, we build a collection of 74 models for human gut bacteria and test their ability to reproduce growth on a set of experimentally defined media. Finally, we create a database of 5587 bacterial models and demonstrate its potential for fast generation of microbial community models. Overall, CarveMe provides an open-source and user-friendly tool towards broadening the use of metabolic modeling in studying microbial species and communities.

Laboratory evolution reveals regulatory and metabolic trade-offs of glycerol utilization in Saccharomyces cerevisiae.

Most microbial species, including model eukaryote Saccharomyces cerevisiae, possess genetic capability to utilize many alternative nutrient sources. Yet, it remains an open question whether these manifest into assimilatory phenotypes. Despite possessing all necessary pathways, S. cerevisiae grows poorly or not at all when glycerol is the sole carbon source. Here we discover, through multiple evolved lineages, genetic determinants underlying glycerol catabolism and the associated fitness trade-offs. Most evolved lineages adapted through mutations in the HOG pathway, but showed hampered osmotolerance. In the other lineages, we find that only three mutations cause the improved phenotype. One of these contributes counter-intuitively by decoupling the TCA cycle from oxidative phosphorylation, and thereby hampers ethanol utilization. Transcriptomics, proteomics and metabolomics analysis of the re-engineered strains affirmed the causality of the three mutations at molecular level. Introduction of these mutations resulted in improved glycerol utilization also in industrial strains. Our findings not only have a direct relevance for improving glycerol-based bioprocesses, but also illustrate how a metabolic pathway can remain unexploited due to fitness trade-offs in other, ecologically important, traits.